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This Letter describes the first determination of bounds on the CP- violation parameter 2f3 s using 
B° decays in which the flavor of the bottom meson at production is identified. The result is based 
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observed data is 15%, corresponding to 1.5 Gaussian standard deviations. 

Dedicated to the memory of our dear friend and colleague, Michael P. Schmidt. 

PACS numbers: 13.25.Hw, 14.40.Nd, 14.65. Fy 



* With visitors from "University of Athens, 15784 Athens, Greece, j, tt . ., , . ., c , „ nor^n w j -j o 

J ' ' ' "Umvcrsidad Autonoma of Madrid, E-28049 Madrid, Spain, 



4 



The accurate determination of charge-conjugation- 
parity (CP) violation in meson systems has been one 
of the goals of particle physics since the effect was first 
discovered in neutral kaon decays in 1964 Stan- 
dard model CP-violating effects are described through 
the Cabibbo-Kobayashi-Maskawa (CKM) mechanism [2j , 
which has proved to be extremely successful in describing 
the phenomenology of CP violation in B° and B + decays 
in the past decade However, comparable experimen- 
tal knowledge of B® decays has been lacking. 

In the B® system, the mass eigenstates B® L and B° sH 
are admixtures of the flavor eigenstates B® and B^. This 
causes oscillations between the B® and B® states with 
a frequency proportional to the mass difference of the 
mass eigenstates, Am s = ran — mf,. In the standard 
model this effect is explained in terms of second-order 
weak processes involving virtual massive particles that 
provide a transition amplitude between the B® and B® 
states. The magnitude of this mixing amplitude is pro- 
portional to the oscillation frequency, while its phase, 
responsible for CP violation in — > J/tp </> decays, is 

-2/9™ = -2arg(-pJ$|) 0, where V i3 are the el- 
ements of the CKM quark mixing matrix. The pres- 
ence of physics beyond the standard model could con- 
tribute additional processes and modify the magnitude 
or the phase of the mixing amplitude. The recent precise 
determination of the oscillation frequency [H indicates 
that contributions of new physics to the magnitude, if 
any, are extremely small [6]. Global fits of experimen- 
tal data tightly constrain the CP phase to small values 
in the context of the standard model, 2/3f M » 0.04 Q. 
However, new physics may contribute significantly larger 
values [y, [j| . The observed CP phase can be expressed 
as 2/3 s = 2/3f M - $ p , where <j)f p is due to the addi- 
tional processes. The decay-width difference between the 
mass eigenstates, Ar = Tl — Th, is also sensitive to the 



same new physics phase. If 



lNP 



> 2/3* M , we expect 



Ar = 2|ri 2 |cos(2/3 s ) [a], where |r 12 1 is the off-diagonal 
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element of the B Q S -B® decay matrix from the Schroedinger 
equation describing the time evolution of B® mesons Q ■ 
Recent studies of B° — > J/ip4> decays without identifi- 
cation of the initial flavor of the B® meson 0, [l(| have 
provided information on Ar and have some limited sen- 
sitivity to the CP phase. 

In this Letter we present the first study of the B® — » 
J/ip<fi decay [111 in which the initial state of the B2 me- 



son (i.e. whether it is produced as B® or its anti-particle 
B®) is identified in a process known as "flavor tagging". 
Such information is necessary to separate the time evo- 
lution of mesons produced as B° or B®. By relating 
this time development with the CP eigenvalue of the fi- 
nal states, which is accessible through the angular dis- 
tributions of the J/ip and <f> mesons, we obtain direct 
sensitivity to the CP-violating phase. This phase enters 
the time-development with terms proportional to both 

cos(2/3 s )| and sin(2/3 s ). Analyses of B° s — > J/ip<f> decays 
that do not use flavor tagging are primarily sensitive to 

cos(2/3 s )| and | sin(2/3 s )|, leading to a four-fold ambigu- 
ity in the determination of 2f3 s [ij |T(| • 

This measurement uses 1.35 fb _1 of data collected by 
the CDF experiment at the Fermilab Tevatron between 
February 2002 and September 2006. The CDF II detector 
is described in detail in Ref. 12j. Detector sub-systems 



relevant for this analysis are described briefly here. The 
tracking system is composed of silicon micro-strip de- 
tectors surrounded by a multi-wire drift chamber. The 
drift chamber provides tracking information and charged 
particle identification through the measurement of spe- 
cific ionization energy loss (dE/dx). A time-of- flight 
(TOF) detector provides additional particle identifica- 
tion. These detectors are immersed within a 1.4 T axial 
magnetic field. Electromagnetic and hadronic calorime- 
ters surround the solenoid. At the outermost radial ex- 
tent of the detector, muons are detected in planes of 
multi-wire drift chambers and scintillators. The data 
used were collected with a di-muon trigger which pref- 
erentially selects events containing J/ip — > de- 
cays [lif . 

We reconstruct the B° — ► J/tj) 4> decay from the decays 
J/tp — > and <fi ^ K + K~ and require these final 

state particles to originate from a common point. We use 
an artificial neural network (ANN) [l3| to separate B® — > 
J/4>4> signal from background. In the ANN training, we 
consider the following variables: particle identification 
of kaons using the TOF and dE/dx, the component of 
momenta of the B® and <f> candidates transverse to the 
proton beam direction, the invariant mass of the <p candi- 
date, and the quality of a kinematic fit to the trajectories 
of the final state particles. We have trained the ANN 
with signal events from simulated data that are passed 
through the standard GEANT-based 14( simulation of the 
CDF II detector [HI and are reconstructed as in real 
data. We use B° — > J/ip4> mass sideband candidates, 
defined as those having m(J/V><£) G [5.1820,5.2142] U 
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[5.3430, 5.3752] GeV/c 2 , as the background sample in the 
ANN training. Applying the selection on the output vari- 
able of the ANN, we observe 2,019 ± 73 B° -> J/ip^ 
signal events with a signal to background ratio of ap- 
proximately one. The invariant J/tpfi mass distribution 
is shown in Fig. [T] An event-specific primary interaction 
point is used in the calculation of the proper decay time, 
t = m(B s )L xy (B°)/p T (B°), where L xy {B° s ) is the dis- 
tance from the primary vertex to the B® — > J/tp (f> decay 
vertex projected onto the momentum of the B® in the 
plane transverse to the proton beam direction, m(P") is 
the mass of the B® meson [jj, and pt(B®) is its measured 
transverse momentum. 




5.45 

Mass(JA|/<|>) (GeV/c 2 ) 



FIG. 1: Invariant p + fx~ K + K~ mass distribution with the 
fit projection overlaid. The vertical lines indicate the mass 
sideband regions. 



The orbital angular momenta of the vector (spin 1) 
mesons, J/tp and <f>, produced in the decay of the pseu- 
doscalar (spin 0) B° B meson, are used to distinguish the 
CP-even S- and D-wave final states from the CP-odd P- 
wave final state. We measure the decay angles 0t, (/>r, 
and ipT j defined in Ref. [§] , in the transversity basis [lj| • 
The transverse linear polarization amplitudes at t = 0, 
Au and A±, correspond to CP even and CP odd final 
states, respectively. The longitudinal polarization am- 
plitude Aq corresponds to a CP even final state. The 
polarization amplitudes are required to satisfy the con- 
dition \A \ 2 + \A l{ \ 2 + \A ± \ 2 = 1. 

In order to separate the time development of the 
B° meson from that of the B® meson, we identify the 
flavor of the B a a or B° s meson at the time of production 
by means of flavor tagging. Two independent types of 
flavor tags are used, each exploiting specific features of 
the production of b quarks at the Tevatron, where they 
are mostly produced as bb pairs. The first type of flavor 
tag infers the production flavor of the B° s or B° s meson 
from the decay products of the b hadron produced by the 
other b quark in the event. This is known as an opposite- 
side flavor tag (OST). The OST decisions are based on 
the charge of muons or electrons from semileptonic B de- 



cays [17|,ll8j or the net charge of the opposite-side jet [191 ]. 
If multiple tags are available for an event, the decision 
from the highest dilution flavor tag is chosen [20]. The 
tag dilution V, defined by the probability to correctly 
tag a candidate P tag = (1 + T>)/2, is estimated for each 
event. The calibration of the OST dilution is determined 
from B+ -> J/ip K+ and B° -> J/tp K*° decays. The 
second type of flavor tag identifies the flavor of the re- 
constructed P° or P° meson at production by correlating 
it with the charge of an associated kaon arising from frag- 
mentation processes 21| , referred to as a same-side kaon 
tag (SSKT). The SSKT algorithm and its dilution cali- 
bration on simulated data are described in Ref. [22J . The 
average dilution is (11 ± 2)% for the OST and (27 ± 4)% 
for the SSKT, where the uncertainties contain both sta- 
tistical and systematic effects. The measured efficiencies 
for a candidate to be tagged are (96 ± 1)% for the OST 
and (50 ± 1)% for the SSKT. 

An unbinned maximum likelihood fit is performed to 
extract the parameters of interest, 2(3 S and AT, plus nui- 
sance parameters to the measurement, which include the 
signal fraction f s , the mean B° s width T = {Tl + r#)/2, 
the mixing frequency Am Sl the magnitudes of the po- 
larization amplitudes |A)| 2 , I^||| 2 j anc * \A ±\ 2 , and the 
strong phases = arg(A|A ) and 5± = arg(^Ao). 
The fit uses information on the reconstructed B® candi- 
date mass m and its uncertainty a m , the B° s candidate 
proper decay time t and its uncertainty a t , the transver- 
sity angles p = {cos 8t, 4>t, cos ipr}, and tag informa- 
tion T> and £, where T> is the event-specific dilution and 
£ = { — 1,0, +1} is the tag decision, in which +1 corre- 
sponds to a candidate tagged as B®, — 1 to a and 
to an untagged candidate. The single-event likelihood 
is described in terms of signal (P s ) and background (P&) 
probability distribution functions (PDFs) as 



f s P s {m\a m )P s {t,p^\V,a t )P s {a t )P s {V) 
+ (1 - f s )P b (m)P b (t\a t )P b (p)P b (a t )P b (V). 



(1) 



The signal mass PDF P s (m\a m ) is parameterized as a sin- 
gle Gaussian with a standard deviation determined sep- 
arately for each candidate, while the background mass 
PDF, Pf,(m), is parameterized as a first order polyno- 
mial. The distributions of the decay time uncertainty 
and the event-specific dilution are observed to be differ- 
ent in signal and background, so we include their PDFs 
explicitly in the likelihood. The signal PDFs P s {<Jt) and 
P s (V) are determined from sideband-subtracted data dis- 
tributions, while the background PDFs P b (at) and P b {T>) 
are determined from the J/ip(f> invariant mass sidebands. 
The PDFs of the decay time uncertainties, P s (<Jt) and 
Ph(cr t ), are described with a sum of Gamma function dis- 
tributions, while the dilution PDFs P s (2?) and P b (T>) are 
included as histograms that have been extracted from 
data. 

The time and angular dependence of the signal PDF 
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P s (t, p, £, \T>, at) for a single flavor tag can be written in 
terms of two PDFs, P for B° and P for _B° as 



P s {t,p,(,\V,a t ) 



2 

1-£D 



P(t,p\a t )e(p) 
P(t,p\a t )e(p), 



(2) 



which is trivially extended in the case of two indepen- 
dent flavor tags (OST and SSKT). The detector accep- 
tance effects on the transversity angle distributions, e(p), 
are modeled with B® — * J/ip<j> simulated data. Three- 
dimensional joint distributions of the transversity angles 
are used to determine e(p), in order to correctly account 
for any dependencies among the angles. The time and 
angular probabilities for B° s can be expressed as 



d 4 P(t,p) 
dtdp 



cx |Ao| 2 T + / 1 (p) + |A||| 2 T + / 2 (p) 



+ lA^T-hifi + lA^A^U+Uip) 

+ |A ||A|||cos(5||)T + / 5 (p) 

+ \Ao\\A ± \V+f 6 (fl, (3) 

where the functions fi(p) ■ . ■ fe(p) are defined in Ref. 
The probability P for B® is obtained by substituting 
U+ — > U- and V+ — > V_. The time-dependent term 
7± is defined as 



T±=e 



-n 



[cosh(Ari/2) T cos(2/3 s ) sinh(AK/2) 
=F 77sin(2/3 s )sin(Am s t)] , 



where 77 = +1 for P and —1 for P. The other time- 
dependent terms are defined as 



U± = ±e 



-rt 



x [sin((5j_ — S\\) cos(Am s t) 

— cos(S± — 6\\ ) cos(2/3 s ) sin(Am s i) 

± cos(<5 ± - 5\\) sin(2/3 s ) sinh(AR/2)] 



V± = ±e x [sin(5x) cos(Am s f) 

— cos(<5^) cos(2/3 s ) sin(Am s i) 
± cos(<5_l) sin(2/3 s ) sinh(AR/2)] . 

These relations assume that there is no direct CP viola- 
tion in the system. The time-dependence is convolved 
with a Gaussian proper time resolution function with 
standard deviation a t , which is adjusted by an overall 
calibration factor determined from the fit using promptly 
decaying background candidates. The average of the res- 
olution function is 0.1 ps, with a root-mean-square devi- 
ation of 0.04 ps. 

We model the lifetime PDF for the background, 
Pb{t\at), with a delta function at t = 0, a single neg- 
ative exponential, and two positive exponentials, all of 
which are convolved with the Gaussian resolution func- 
tion. The background angular PDFs are factorized, 
Pb(p) = Pb(cos9T)Pb(<PT)Pb(cosipT), and are obtained 
using B® mass sidebands events. 



Possible asymmetries between the tagging rate and 
dilution of B® and B® mesons have been studied with 
control samples and found to be statistically insignifi- 
cant. We allow important sources of systematic uncer- 
tainty, such as the determination of overall calibration 
factors associated with the proper decay time resolu- 
tion and the dilutions, to float in the fit. The mixing 
frequency Am s — 17.77 ± 0.12 ps^ 1 is constrained in 
the fit within the experimental uncertainties 0. Sys- 
tematic uncertainties coming from alignment, detector 
sculpting, background angular distributions, decays from 
other B mesons, the modeling of signal and background 
are found to have a negligible effect on the determination 
of both Ar and f3 s relative to statistical uncertainties. 

An exact symmetry is present in the signal probability 
distribution, as can be seen in Eq. ([3]), which is invariant 
under the simultaneous transformation (2/3 s — > ir — 2p s , 
Ar — > — Ar, Si\ — > 2tt—S», and(5_L — > tt— S±). This causes 
the likelihood function to have two minima. This sym- 
metry can be removed by restricting any of the above pa- 
rameters within appropriate ranges. However, even after 
removal of the exact symmetry, approximate symmetries 
remain, producing local minima. Since the log-likelihood 
function is non-parabolic, we cannot meaningfully quote 
point estimates. Instead we choose to construct a confi- 
dence region in the 2j3 s — Ar plane. 

We use the Feldman-Cousins likelihood ratio order- 
ing [23[ to determine the confidence level (CL) for a 
20 x 40 grid evenly spaced in 2/3 s <G [— 7r/2, 3tt/2] and 
Ar G [—0.7,0.7]. The other parameters in the fit are 
treated as nuisance parameters (e.g. B® mean width, 
transversity amplitudes, strong phases) [24|. The cover- 
age against deviations of the nuisance parameters from 
the measured values is confirmed by randomly sampling 
the nuisance parameter space within ±5ct of the values 
determined from the fit to data. The 68% and 95% confi- 
dence regions obtained are shown in Fig. [2] The solution 
centered inO < 2/3 s < ir/2 and Ar > corresponds 
to cos(<5_l) < and cos(<5j_ — S\\) > 0, while the oppo- 
site is true for the solution centered in it/2 < f3 s < ir 
and Ar < 0. Assuming the standard model predicted 
values of 2/3 s = 0.04 and Ar = 0.096 ps" 1 [8], the prob- 
ability to observe a likelihood ratio equal to or higher 
than what is observed in data is 15%. Additionally, 
we present a Feldman-Cousins confidence interval of 2(3 Sl 
where Ar is treated as a nuisance parameter, and find 
that 2/3 s € [0.32, 2.82] at the 68% confidence level. The 
CP phase 2(3 S , AT, T, and the linear polarization am- 
plitudes are consistent with those measured in Ref. 0]. 

We also exploit current experimental and theoreti- 
cal information to extract tighter bounds on the CP- 
violating phase. By applying the constraint | Ti2 1 = 
0.048 ± 0.018 [8] in the relation Ar = 2|r 12 | cos(2/3 s ), 
we obtain 2(3 S G [0.24, 1.36] U [1.78,2.90] at the 68% CL. 
If we additionally constrain the strong phases Su and S± 
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FIG. 2: Feldman-Cousins confidence region in the 2f3 s — Ar 
plane, where the standard model favored point is shown with 
error bars 0. The intersection of the horizontal and vertical 
dotted lines indicates the reflection symmetry in the 2f3 a — Ar 
plane. 

to the results from B° — > J/ijjK*° decays [25| and the 
B® mean width to the world average B h width [|, we 
find 2(3 S 6 [0.40, 1.20] at the 68% CL. 

In summary we present confidence bounds on the CP- 
violation parameter 2f3 s and the width difference Ar 
from the first measurement of B® — > J/ip 4> decays using 
flavor tagging. Assuming the standard model predicted 
values of 2(3 S = 0.04 and Ar = 0.096 ps -1 , the proba- 
bility of a deviation as large as the level of the observed 
data is 15%, which corresponds to 1.5 Gaussian standard 
deviations. Treating AT instead as a nuisance parameter 
and fitting only for 2f3 s , we find that 2/3 s 6 [0.32,2.82] 
at the 68% confidence level. The presented experimental 
bounds restrict the knowledge of 2(3 S to two of the four 
solutions allowed in measurements that do not use flavor 
tagging jjj [l(| and improve the overall knowledge of this 
parameter. 
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